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(-H , We obtain novel nonlinear Schrodinger-Pauli equations through a formal non-relativistic 

f^ I limit of appropriately constructed nonlinear Dirac equations. This procedure automat- 

ically provides a physical regularisation of potential singularities brought forward by 
the nonlinear terms and suggests how to regularise previous equations studied in the 
literature. The enhancement of contributions coming from the regularised singularities 
suggests that the obtained equations might be useful for future precision tests of quantum 
nonlinearity. 
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1. Introduction 



Several nonlinear extensions of Schrodinger's equation have been constructed to 
probe the accuracy of quantum linearity. ^"^ For example, Weinberg proposed a 
class of equations which were then used in several experimental tests, see^"^ and 

QQ \ references therein. The results indicated that any potential non-linearity in those 

^3 ' systems had to be smaller than some bounds. 

Ignoring external fields, the nonlinear Schrodinger equations may be written in 
the form 



S' ih^^^ = -^v^^ + fNRm^ (1) 

where the nonlinearity Jnr depends in general on the wavefunction, its conjugate 
and their derivatives. Jnr may be written as a ratio of two terms, Nijp)/ Diij)), with 
equal factors of ip in the numerator and denominator to keep the scale invariance 
"0 — >■ Ai/), A a constant, of the linear Schrodinger equation. This scale invariance 
ensures the wavefunction can be freely normalised. The denominator is typically a 
monomial in ?/;*?/' so that the nonlinear term may be made separable for indepen- 
dent systems. As the nonlinearity must be weak on phenomenological grounds, the 
solutions of the linear equation must be very close to some solutions of the modified 
equation. However, any solutions of the linear equation that have nodes would make 
-D(^) vanish at some points. Thus, the nonlinearity would generally be singular and 
ill-defined at those points. 
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Weinberg,'^ discusses classes of nonlinear Schrodinger-Pauli equations where the 
nonlinearity turns out to be finite at the nodes because the numerator vanishes 
faster than the denominator there. However this will not happen for general classes 
of nonlinearities where N{ip) has derivatives, such as the equations studied in Ref.^ 

It was suggested in Ref.^ that quantum nonlinearity might be linked to the 
breaking of space-time symmetry. This idea was supported by a study in Ref.^° in 
the relativistic regime: That is, a deviation from quantum linearity is associated 
with a violation of Lorentz symmetry. Thus, in this paper, we discuss how to con- 
struct novel classes of nonlinear Schrodinger-Pauli equations, which have the above- 
mentioned scale invariance, starting from nonlinear Dirac equations. However, we 
will keep our option open by considering both Lorentz invariant and Lorentz vio- 
lating nonlinear Dirac equations. ^^ 

As we shall see, our procedure of obtaining the Scrodinger-Pauli equations by 
taking the limit of relativistic equations has the advantage of indicating a natural 
and physical regularisation of the singularities. We remark that we focus on genuine 
nonlinear Dirac equations that cannot be linearised by performing a nonlinear gauge 
transformation. ^ ^ 

In the next section, we discuss, in general terms, the formal non-relativistic limit 
of a subset of nonlinear Dirac equations constructed in Ref.-'^"'^ For conciseness, in this 
paper we only consider the case where F = fl in ^, I being the identity matrix 
in spinor space. Explicit examples of the lowest order nonlinearities, corresponding 
to one factor of ^p*ip in D{ip) are exhibited in Section [31 other cases being similarly 
handled. The singularity resolution is discussed in Section 2] and we end with a 
discussion in Section [5] 

We note in passing that nonlinear Schrodinger equations of other types have 
been constructed from Levy-Leblond's "non-relativistic Dirac equation" which is 
itself the non-relativistic limit of the usual Dirac equation. ^^'^^ 

2. Non-Relativistic Limit 

We start from nonlinear Dirac equations of the form 

(ift7^a^ ~mc + eF)iP^O, (2) 

where F — F[ip, tp) — fl and where we have made the small parameter e explicit. 
We demand that F has certain properties so that desirable characteristics of the 
linear Dirac equation, such as locality, conservation of probability, separability and 
invariance under tp — >■ Xip, are retained (we are adopting the standard kinemat- 
ical structure of quantum mechanics, in particular the standard inner product). 
The other symbols in ([2]), such as those for the gamma matrices, have their usual 
meanings; our conventions are similar to those in the textbook^* and in Rcf.^^ 
In Hamiltonian form the equation is 

ih—'iP = {ihca ■ V + I3mc^ - ecfif) ip (3) 
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where a' = j^j'^ and /3 = 7°. It maybe be decomposed into two equations by 
introducing upper and lower components of the wavefunction, 

^^ /'¥''\e-^mc^tA (4) 

where the rest energy has been extracted as it is the largest component in the 
non-relativistic limit. We adopt the standard textbook procedure in obtaining the 
leading non-relativistic limit, but for clarity we repeat some steps below. In order 
to make the algebra manageable, we simply take 1/c to be the same order of mag- 
nitude as the nonlinearity scale e and keep only the leading nonlinear term in the 
standard non-relativistic expansion. Thus we can isolate the leading order nonlin- 
ear contribution. However, in realistic applications, e will be much smaller than 1/c: 
This will introduce higher order, 1/c, terms which do not affect the leading order 
nonlinear contribution. 

Substituting ^ into ^ we get 

(As in the usual textbook procedure, the ansatz (|4]) removes the mass term for the 
upper component). From the lower component of ([5]) we have, 

^ iha- ■ V(p ih dx efx /gx 

2mc 2mc'^ dt 2mc 

Let Xo = ''^f^e "^ ■ Then expanding © about xo, we obtain ^ = Xo + O (^, ■^). 
That is, X is the same as that in the linear theory. Substituting ([6]) into the upper 
component of ([5]), we arrive at 

'^dt^^-L"^'^-''^'"'^ ^^) 

where f^R means that the state dependence of / has been simplified using ((H [S]) 
and higher order terms dropped. Below we look at some explicit examples. 

3. Examples 

3.1. Lorentz invariant f with one derivative 

A Lorentz invariant / with one derivative and which is odd under the parity trans- 
formation is 

/i = ^%f , (8) 

where jg = V'7^75V' is the usual chiral current. The non-relativistic limit is 

dip h^V'^(p V • iip^cTip) 

at 2ni \Lp\ 

The factor V • {^Lp^ crip] appears often in parity odd equations;^^ it couples the spin 
components of the two-component spinor. 
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3.2. Lorentz invariant f with two derivatives 

For an example of a Lorentz invariant / with two derivatives consider 

■h -r-f , (J-UJ 

ipy; 

where e and S are two independent small parameters taken to be of order 1/c below. 
The non-relativistic equation is 

.,d(p fi?\7'^ip / mcS\ If ( Simc j.dtf ( df^ 

~dt ~ 2^ [^ + '2f^) + M^ t fi~ [^ ~dt~ [~dr 

+ ec(vVV)+'^c(V(^t) .(v^)} , (11) 



3.3. Lorentz violating, parity even f 

Lorentz violating non-linear Dirac equations are of some interest. ^'^"'^'"'^^^^^ An ex- 
ample of such an / with no derivatives and even under parity is 

where A^^ is a constant vector background field. The non-relativistic limit is 

..dip ft2y2 zhipA-[if^Vif-{Wf^)if] 

*^li7 o cAoip + - ^ —r^ '-^. (13) 

ot 2m 2m |(/?p 

3.4. Lorentz violating, parity odd f 

A Lorentz violating / which is odd under parity is 

U-A,tfpt, (14) 

The non-relativistic equation is 

dip ffV'^ip cip'^A-cnp Aoihif [ip^cr -Vifi — {\7ip^) ■ (Tip] 



4. Apparent Singularities 

From the above examples, we see the appearance of the following structures in the 
non-linear Schrodinger-Pauli equations, 

^ — — i — f^ — I '^ — i — n ' ^ ^ — i — 12 — ■ y'-^' 

Clearly, at the nodes of ip, these forms are singular. However, we can avoid these 
singularities in a natural way. For our nonlinear Dirac equations, ^^ the nonlinear- 
ities have the common structure /j^'Z, , the n — 1 case being discussed here. In 
terms of the two component spinors this is i i^^i i ^ where the lower (small) compo- 
nent contribution |xp is usually dropped in the non-relativistic limit. However, at 
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the nodes of </?, we must keep the small component in the denominator. This regu- 
lates the above mentioned singularity for the following reason: From ([B|), the lower 
component is proportional to the slope of Lp (i.e. V(/?), which is unlikely to vanish 
simultaneously at the nodes except for special cases. In such extreme cases, one 
would need to retain the smaller terms (higher order in 1/c) in the non-relativistic 
expansion of the denominator. 

For the specific examples illustrated above, the replacement |(^|^ — > \lp\'^ — |xP 
in the denominator makes X = Z = at a, node of (p while Y becomes finite and 
actually enhanced because of the small denominator. Note that, at the level of the 
equation of motion, there is an extra factor of ip which multiplies the nonlinearity 
/. It is clear that X and Z contributions in the equation of motion are not singular 
even at the nodes but the Y contribution is, unless regularised as discussed above. 

So far we have discussed singularities in / and at the level of equations of motion. 
As for observables, let us consider shifts in the energy levels given by first-order 
perturbation theory, 

6E= f d^x <tp\F\ip>^ f d^x\ip\^f{ip) (17) 

where the unperturbed (linear equation) wavefunctions are used. We see that the 
X, Y, Z structures give finite shifts. Singularities will appear in n > 2 classes of 
nonlinearities discussed in Ref,^^ two examples of which are given by 

|(^|2|(^|2 
[(V^t) . (V^)] (y.tv2^) 

It is clear that the energy shifts will be singular for such terms unless the regulari- 
sation is implemented. 

The above discussion has ignored external potentials which must be included in 
realistic experiments. For example, in the presence of an external gauge field and for 

a particular spin component ip = { „ ] fo, the lower component is modified from 



W^YZ^ ^' ^ ' I.I2M2 -■ (19) 



,0 
its previous form xo to become 



Xo^T^ij. A. ) <^o ■ (20) 

2mc \oz c ' 



Setting ipo = g[x — Xq) near a node we have 

,2_ ^' 



IXoI 



4y^2g2 



''^^ 2'-A.g^+'4Al9' 



dz I c ' dz c? 



(21) 



In this case, at the node of ip, |xoP has exactly the same form as the case when the 
gauge field is absent. 

We can see that the contributions from nonlinear effects are largest (if non-zero) 
at the nodes. This suggests that future tests for quantum nonlinearity should focus 
on systems containing nodes in their wavefunctions. 
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We remark that the nonUnear equations discussed in Ref.® have been appUed to 
the hydrogen system^* but the physical consequences of singularities at the nodes 
of wavefunctions {ip -^ 0) was not discussed. 

5. Discussion 

We have illustrated how to obtain novel classes of nonlinear Schrodinger-Pauli equa- 
tions starting from the nonlinear Dirac equations constructed in Ref,^^ the latter 
equations themselves being more general than previous constructions. ^^"^^ For ex- 
ample, we have cases where the time-derivatives appear in the nonlinearity, and cases 
where the two components of the spinor are coupled through parity violation. We 
remark that probability is conserved for all of our non-relativistic equations. Also, 
the equations that are descended from Lorentz covariant equations are Galilean 
invariant. 

An interesting point to note is that certain Lorentz-violating nonlinear Dirac 
equations have non-relativistic limits that are Galilean invariant. For example, for 
/3, if the background field has only a time component, the leading non-relativistic 
limit actually becomes linear and invariant under Galilean transformations. For f^, 
choosing a space-like background field will cause the non-relativistic equation to be 
still nonlinear but invariant under Galilean transformations. 

We had taken the nonlinearity parameter e to be the same order of magnitude as 
1/c for ease of power counting, as our main aim was to isolate the leading nonlinear 
structure in the formal non-relativistic limit. We saw that potential singularities 
in the Schrodinger-Pauli equations are regularised by keeping the sublcading lower 
components of the four component Dirac spinor in the denominators of the non- 
linear terms. Thus, physically, it is the relativistic corrections that regulate the 
singularities. Precisely at a node, if the numerator is is nonzero, the nonlinearity is 
actually enhanced by the small denominator. 

The situation here is qualitatively similar to a previous study^^ of an 
information-theoretic motivated nonlinear Schrodinger equation,^ where the con- 
tribution to energy shifts from states with nodes was enhanced relative to states 
which had no nodes. Note also that in replacing the potentially singular denom- 
inator |(^p by \(p\'^ — Ixl^ as in Section SI one has introduced an infinite number 
of derivatives, through a formal expansion of the denominator, into the nonlinear 
terms even though we had started with a finite number of derivatives. This again is 
qualitatively similar to the situation with the information-theoretic nonlinearity^ 

In actual applications, such as tests of quantum linearity, one would have to set 
e much smaller than 1/c in the constructed nonlinear Schrodinger-Pauli equations 
even though they were formally derived from the nonlinear Dirac equations assuming 
e - 1/c. 

The main suggestion from this study is that future precision low-energy exper- 
iments, probing deviations from quantum linearity, should focus on systems which 
have nodes in their limiting linear wavefunctions. It is there that the nonlinearity, 
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if nonzero, will be enhanced. 
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